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' Abstract 

The properties of measuring systems of the Extensive Air Showers arrival directions (EAS 
goniometers) , consisting of some sets of widely separated detectors registering the moments of 
EAS frontal passage, are considered. The advantage of volumetric goniometers in comparison 
with the commonly used flat ones is shown. The proper points selection method for detectors 
spatial arranging is suggested, providing the best accessible accuracy of the EAS arrival 
direction estimation within the given detectors' number and installation overall dimensions. 
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1 Introduction 



Since the discovery of cosmic rays (CR) the problem of large-scale correlations has been almost 

. disregarded. This is because charged particles, composing a main part of CR, must be deflected 

in the galactic magnetic field throughout their long journeys to the Earth, resulting in random 

, injections into atmosphere of CR particle pairs, coinciding both in time and in direction. However, 

' there exist some processes, engendering such pairs and groups of correlated primary particles in 

OO ' 

' cosmic radiation. There may be, e.g., groups from the 7-bursts, the product 7-rays from super 

high energy collisions in the interstellar substance (especially in the immediate vicinity of CR 
generation areas), or the products of CR nuclei disintegration in the solar photon field - i.e. the 
effect presented by Gerasimova and Zatsepin [Hll], etc. 

All this effects are very rare, but during the last years there appeared some observational 
evidence of their existence [31 HI • All of them were performed by observation of correlated 
Extensive Air Showers (EAS). Actually, only both the registration times' coincidence and the 
coincidence of two EAS cores' directions with account of spatial positions of registration points 
on the Earth's surface, can state as a partial warranty of their genetic relationship 

The problem can be solved in principle by the pool of modern complicated installations, 
investigating the 7-initiated EAS, such as the MAGIC in the Canary Islands [6] or VERITAS 
in southern Arizona [7J. These installations detect and log out the showers by observation of 
Cherenkov light generated in the atmosphere by the quick charged particles composing the EAS. 
The arrival direction measurement precision of 0.15° is expected for individual 7-quanta. 

However, these extremely perfect installations are rarely disposed on the Globe. The active 
observation time is limited by the atmosphere state condition. The viewing angle usually is 
moderate, e.g. the VERITAS array ^ observes concurrently the sky area of 3.5° only. Thus it is 
not expectable for these expensive installations to constitute some effective pool for uninterrupted 
monitoring of rare events of widely spaced genetically related showers. 



It is safe to say that the widespread network of small traditional ground-based installations, 
recording both the EAS event absolute astronomical time (UT) and the shower arrival direction 
on the local upper hemisphere (and measuring the EAS energy too, if possible), would make up 
more suitable tool for the problem. The installations of this type, being comparatively simple and 
inexpensive, would form reasonable equipment for university teems and research groups round 
the Globe, constituting the synchronized united network. 

Now, as a matter of notational convenience, let us term any ground-based installation, 
intended for the EAS arrival direction measurement, as EAS goniometer. 

The LAAS group [9j has performed the network observation of EAS with energy 
(10^^ -ir 10^^)eV, using ten independent scintillator goniometers arranged over a very large area. 
Especially this last group has investigated the EAS pairs coincidence not only by the times reg- 
istered (using the Global Positioning System GPS) at the network stations, but also by the EAS 
arrival directions coincidence. This group has declared the attained measurement accuracy of 
EAS arrival direction to be about 10° for the zenith angles less then 45°. This attainment is too 
small for the reliable paternal affiliation in the pairs observed. The analysis of the typical station 
construction of the network [9j gives the opportunity to ascertain the possibility of much greater 
angular accuracy achievement even with the same number of detectors as in the EAS goniometer 
used there. It becomes possible to get an acceptable accuracy even for big zenith angles. This 
accuracy growth is very important for any network of remote stations as far as the potentially 
correlated showers prove to be close to the horizon for big enough angular distances between the 
stations. Just for the last cases the extensively used flat EAS goniometers with poor accuracy at 
large zenith angles are especially objectionable. 

2 Volumetric EAS goniometer. Common case. 

Generally the set of any Nq ionizing radiation detectors, arbitrarily distributed in the 3-space, 
can be used as the volumetric EAS goniometer. Some "triggering structure" for EAS discovering 
is implied. This structure has to send a trigger signal to the measuring part of installation to 
start a time reckoning of the pulses from EAS goniometer detectors. In particular, the EAS 
goniometer detectors themselves can be used for this purpose. The signals from all detectors 
have to be delayed for the common period, t^eh with respect to the real moments of the pulses 
origins. This delay period has to be such that at any possible case of the shower arrival direction 
all A'o detector signals would hit the measuring part of installation later then the trigger signal. 
The measuring part itself records the differences of signal arrival times with respect to the trigger 
signal hit time as well as the absolute astronomic time (UT) of trigger signal obtained by means 
of GPS. Only this set of recorded numbers has to be used later (off line) to estimate the direction 
and the shower front passing time through the installation's coordinate system origin. 

Hereinafter all delay periods are measured in distance unites , that is to say their clock esti- 
mations are multiplied in advance by the EAS front velocity, which is taken to be approximately 
equal to the light velocity. 

Let us fix the rectangular coordinate system for the volumetric EAS goniometer. All 
coordinates, times, their differences and distances are referred hereinafter to the mentioned co- 
ordinate system and are measured in meters. Furthermore, let us designate the coordinates of 
used detector centers as 

ri{x,y,z) i = 1,2, 3... No 



The EAS front plane equation [TU] at arbitrary moment is 



r • n — p = (1) 

Here r is an arbitrary point in the front plane; vector n is the main unit ort of the plane, 
= 1, while the p parameter measures the distance between the plane and coordinate system 
origin. Our aim is to estimate the ort components by the measured times of the front passage 
through the detectors. 

The distances from the detectors to the front plane in any position, specified by p parameter 
value, is determined by the linear relationship 

<5i = r-n-p i = 1,2,3... A^o (2) 

Let us select from the whole family of planes, corresponding to the different moments of 
the shower propagation, the unique plane containing the coordinate origin. The last restriction 
implies the selection of the special plane from the planes' family ([T]), singled out by the condition 
p = . Just this plane will be referred as "the EAS front plane" . 

For this unique EAS front plane the set ([2]) of detectors' distances from the front plane are 

6i = rn (3) 

There can be some negative distances between them in common case, as any of the detector 
points Tj can be disposed on any side the EAS front plane. That is why we shall add an artificial 
and big enough common delay period tdei to both sides of equation ([3]) (for instance, this delay 
can be attained by use of coupling cables with identical and big enough lengths): 

r • n + tdel = 6i + tdel = ttr + U (4) 

Here ti values are the measured by the EAS goniometer installation delay periods of de- 
tectors' signals with respect to trigger signal (all of them are the positive values), while ttr is the 
(unknown) triggering time. 

The decision variables in equation Q are: the common difference of delay period of the 
moment when the EAS front plane passes the coordinate system origin of the installation with 
the triggering time, ts = tdei — Ur, and three nondimensional components of main ort n, i.e. 
the directional cosines of EAS core. The common difference period ts permits to obtain the UT 
of EAS front plane pass through the coordinate system origin (with use of the GPS-measured 
triggering time and common signal delay period known). Since there are 4 decision variables, 
so the requirement on the detectors number Nq > 4 is the solvability condition of equation 
It is possible to get the solution even for three detectors, but it proves to be the singular case 
of system ([3]), demanding a special solution method. This case of "flat" EAS goniometer is 
considered later. Hereon the common case is studied. 

For the case of detector number Nq > 4 the EAS goniometer equations system 

r ■ n + ts = ti (5) 

is overfilled. The Least-Squares Method (LSM) [HI [13] is to be used to determine the solution. 

It is convenient to define new vector ly with 4 rows for the generalized main ort of the 
shower front plane, with the nondimensional 4*^^ component ut to measure the value of desired 
common delay difference. The coordinates of all used detectors constitute the 4 x A'^o matrix p 



with dummy 4*^ column, containing identical 4*^^ component value A = Im for all detectors. The 
measured delay periods constitute the A'^o-row vector r. Now the equation ([5|) takes the form of 
common matrix one: 

pv = T (6) 

As usual [T3] for the overfilled equation set in the frames of LSM, for delay periods 
distributed with dispersion matrix S, let us apply the left multiplication by {p^T,~^) matrix 
to the equation This results in normal equations set for the 4 decision components of the 
generalized direction ort of the EAS front plane 

Au = d (7) 

with the symmetric 4x4 square matrix 

yl = /S-V (8) 

and new right-hand member of equation 

d = p'^^-^T (9) 

It has the unique solution for the nonsingular matrix A. Just due to ([8]) it becomes necessary 
to arrange the detectors set in full 3D space, not in any plane in it. The explicit solution is 

u = A-^d (10) 

or, directly expressed by the delay vector r 

u = Gt (11) 

Here G is the A'^o ^ 4 matrix 

a = A-i/s-i = (/s-V)-V^s-i (12) 



The estimation (jllh of the generalized main ort evaluated here is consistent, unbiased and 
asymptotically normal estimator (i.e. the expectation value of the estimator coincides with the 
true value of the ort), as it is received from the initial data by means of the linear Least Squares 
Method [IS] . 

Essentially, expression (llip is the whole solution in the common case, but the error analysis 
for this solution gives one the possibility to recognize some additional location requirements for 
the final achievement of desired estimation quality of EAS direction. 

At this stage let us assume the coordinates of detectors to be taken as exactly prescribed. 
Practically it means that the detectors' location errors have to be less then the errors of delay 
periods at least by an order of magnitude. 

The generalized direction ort, ly, is linearly connected with the measured delay periods 
by pT]) . So the dispersion matrix, D (of the direction ort v), is connected with the delay disper- 
sion matrix, S, by the linear transformation [1^, too: 

D = G^G^ (13) 



At this stage let us assume as a hypothesis that all delay periods, ti, are identically dis- 
tributed independent quantities with the same dispersion values, cj^. This assumption is very 



close to reality, indeed, as the processes in one delay device do not affect the properties of another 
one. Hence, there are no correlations in the dispersion matrix T, of the delay vector r. On the 
other hand, all delay errors in every signal path originate from the similar reasons, so they can 
be considered as equal on average. Therefore we can accept the relation 

E = cj^I (14) 

Here I is the 4x4 unity matrix. 

In this special case the dispersion matrix of the required vector can be expressed in 
explicit form via the A matrix 

D = A-^ (15) 

So, all volumetric EAS goniometer properties are determined through the detectors arrang- 
ing matrix p and delay dispersion value (see (f8l) (fT4|l ). 

Broadly speaking, it is commonly desirable for the sought quantities (i.e. the generalized 
direction ort v components in our case) to be statistically uncorrelated estimators at least. It is 
desired for the component dispersions to be equal, too. These requirements are strong enough. 
They will give us the possibility to define more accurately the best scheme of detectors arranging 
in the space. 



3 Multitier EAS goniometer 

If it is planed to place the EAS goniometer on the flat horizontal plain, the local vertical 
line becomes the preferential direction in the space and it is natural to orientate the installation 
upon the last one. It is clear that the error isotropy of EAS direction estimations for the upper 
hemisphere is desired. 

Let us consider a multitier scheme of volumetric EAS goniometer. It consists of Nq 
detectors arranged on several horizontal plane levels (tiers). Such a construction prevents the 
possible singularity of equation ([7]). It is clear that the desirable azimuth symmetry of errors 
results in axial symmetry of detectors position on every tier. In the simplest case ( see Fig(T]) 
they are placed uniformly along the circumferences with centers based on a common vertical axis. 

All detectors hereon belong to the goniometer subsystem itself, the trigger subsystem 
is not considered. 

Let us use the rectangular frame of reference with horizontal XOY plane and with OZ axes 
directed upwards along the local vertical line. Axes OX is a polar one for the azimuth angles. 

All detectors are situated on M horizontal levels; a = 1, 2, 3 ... M are indexes of these 

levels. 

Every o-level contains Na detectors. The total number of detectors is 

M 

No = J2Na (16) 

a=l 

The radii of the detectors' positions regarding to the OZ axis are Ra at any a-level, while 
Ha are the levels' heights above the XOY plane. 

At any a-level a separate detectors numeration is established: k = 1,2,3,... Na are the 
indexes of detectors at each level. 



Detectors are represented by points in the 
polygon apexes at any level. The detector 
numbers for the example shown are: 




Figure 1: The detectors' position scheme for the multitier EAS goniometer. 



These detectors divide miiformly corresponding circle of Ra radius with the angular step 

- = l <'^' 

The phase shifts of detectors situated along the a-level circumferences are denoted as ipa- 
Every delay period ta^k of detector with item number k, disposed on the a-level, has to be recorded 
as a result of installation triggering due to EAS front passage through the EAS goniometer. 

It is clear that the a-level detectors are situated in the points with the coordinates: 

{Xa,k = Ra COs{kaa + ^fa), Ya,k = Ra Sm.{kaa + </?a), Ha, A) (18) 

— just these sets of coordinates constitute the A'^o— row matrix p of detectors' positions. 

The circles division mode used above provides the definitive calculation [14J of the A matrix (| 
in the explicit form. 

It is convenient to define the "tier averaging" operation for any set of values Qa representing 
some property of every level. We shall designate this operation by broken brackets: 



< Q > 



def 



1 



M 



" a=l 



(19) 



The explicit form of matrix A turns out to be 

/ < fl2 > 





<R^> 

2 








<H^> X<H> 



(20) 



The matrix is independent of the phase shifts on every level [14j. It proves to be a block- 
diagonal one with only nondiagonal elements corresponding to the 3'^'^ (zenith) and 4*^^ (temporal) 
components of the solution z^. The dispersion matrix D = A^^ has the same structure. 

It turns out than, that the correlation between the zenith and temporal components of the 
generalized ort estimation can be cancelled out by a simple selection of levels heights. These 
heights have to satisfy a relation: 

, M 

<H> = —Y^NaHa = (21) 

Hence, some of them have to be negative, i.e. the origin of installation coordinate system 
must be located above some of the lower tiers. It requires only a simple vertical shift of the initial 
coordinate system. As long as the last operation is a mathematical one and does not require any 
technical change in the installation, it can be accomplished at any circumstances. After this shift 
the matrix A reduces into a diagonal one with new value of < >. It means that the columns 
in the detectors' positions matrix p become orthogonal. 

Hereon we shall suppose that the condition (1211) is fulfilled — it doesn't cost anything! It 
is easy now to present an explicit solution of normal LSM equation ([7]), i.e. the estimation of the 
front plane generalized ort: 

= 2<RC>/<R^> 
vy = 2<RS>/<R^> 

= <HT >/< H^> ^ ^ 

vt = <T>/ \ 

Here are used some "tiered averaged" values, evaluated both from measured delay periods 
of detector signals and the angular coordinates of the detectors: 

k = l 

= j-j- sm{kaa + iPa) (23) 

k=l 

-^a ~ "aT ^ ^"■''^ 
k = l 

The common delay periods' difference tg = t^ei — ttr now can be calculated explicitly 
through the definition: 

ts = X-ut = <T> (24) 

i.e. it is simply an arithmetical mean of all delay periods of all detector signals in the EAS 
goniometer. 

The estimation of ort v dispersion matrix D is shown above for identically distributed 
independent delay periods of detectors' signals (jlSp with the same dispersion values cr^. In the 
case of "orthogonal" multitier installation (i.e. with diagonal A matrix) it reads: 

2^2 a2 ^2 
= A, = n — ; = — ; Dt = ^; (25 



It is strongly desirable to achieve an isotropy of estimation accuracy for all 3D-ort n 
components, too. This aim can be reached if the radii and the heights of all tiers would be 



fitted to satisfy tlie relation < E? >= 2 < >, following the D matrix explicit view (|25|) for 
"orthogonal" goniometers. Hence, the values of tiers' heights would to be of the same order of 
magnitude as the radii used, though it may prove to be difficult for realization. 

According with LSM deductions, the true estimation of dispersion o"^ of delay periods is 
the (corrected) average of squared residual differences: 



M Na 

= TlTr J\ Yl Yl i^a,kVx + Ya^k^y + Za^k^z + Aj^t - ta,fc)^ (26) 

(iVo - 4) 

Here the estimation of the front plane generalized ort (j22p ()23p , the detectors coordinates (|18p 
and measured delay periods must be substituted. Thus the dispersion of measured direction can 
be estimated for every EAS event, but the goniometer must contain strictly more than 4 detectors 
(see (j26]l ). 

The dispersion matrix D = displays (|20 p(j25p some special features of the orthogonal 
EAS goniometer scheme under consideration: 

a) all possible correlations can be eliminated by a simple coordinate shift; 

b) the dispersion values of "horizontal" components of the ort estimated are equal and 
proportionate to l/(A'^o < >); the "vertical" one proportionate to l/(A'^o < >); 

c) the dispersion value of difference of common delay periods does not depend on the 
installation overall dimensions and proportionate to l/(A'^o). 

While the 3Z? — ort n is estimated by means of foregoing procedure (the additional "tem- 
poral" component will be out of consideration hereon!), the problem arises of the corresponding 
spherical angles estimating for the EAS arrival direction, i.e. of azimuth angle (j) and of zenith 
angle 9. 

The direction ort n components are defined through these angles with standard relations: 

Ux = SlliO COS(/> 

Uy = siiiO sincj) (27) 
Uz = cos 6 

It is obvious that n vector ought to be of unit length. This condition may be violated 
exceptionally by the errors in the ort components estimations. Thus it is reasonable and handy 
to prefer the angles calculation method exploiting only the components ratios, as it excludes 
the influence of accidental length variation. So, we accept for computations the special form of 
solution of (j271) : 

cot (/) = rix/riy 



The last are nonlinear expressions, so the usual corrections |13) . depending on the estima- 
tions of the dispersion matrix D components (I25p . must be used for angles estimations. 

The dispersion matrix of the angles obtained is a function of the ort n components' disper- 
sions, derived above ([25|) . Let us calculate the matrix of first derivatives of the angles (p8]) upon 
the ort components {nx,ny,nz)- 



Following the method of error propagation [13] , the dispersion matrix of spherical angles is 

/ Cm 



2 



F^{(l},9)D^F^{^,e) (30) 



Here is the spatial part of the dispersion matrix (j25p for ort v components' estimations. 
For the EAS goniometer with the axial symmetry it results in relations 

^^ = ^%; 'yl = Dh cos" e + D.sin^ e- C^e = 0. (31) 
^ sm 

Here Dh = = Dy. 

In this axially symmetric case the {</>0} — covariation vanishes, both dispersions depend on 
the true zenith angle only. The fast increase of dispersion of azimuth angle estimation is a direct 
sequence of spherical coordinate system singularity: at the — limit the azimuth angle value 
is fundamentally indefinite. The zenith angle value is limited at any case. 



4 Flat EAS goniometer 

Let us investigate the possibilities of flat goniometers by the method used. The only difference 
with 3D case consists in mutual equality of the 3^*^ coordinate of every detector: Za^k = const. 
Hence the equation system ([7D becomes singular; it contains no information about the 3^"^ com- 
ponent of the front plane ort; the equations have no complete solution. 

The way out of the difficulty consists in rejection of the 3'''^ component mentioning from 
the LSM equation system. Let us evaluate now the "horizontal" and "temporal" components 
only, erasing both column and row out of the equation system ([7]). So we can estimate the 
"horizontal" and "temporal" ort components only, being consistent, unbiased and asymptotically 
normal estimator, just as in common case (jlip . This is sufficient for the azimuth angle estimation. 
The value of the 3^"^ component can be reconstructed by a formal way from the unity condition 
for the 3D-ort length: 

"^(n^, Uy) = - {nl + nl) (32) 

In this case the 3^^^ component is not the independently estimated value, but a nonlinear 
function of the estimations of the "horizontal" components. Except the last moment the solution 
of problem is similar to that obtained above for the common case. 

But expressions ()22p and (j32p only nominally solve the problem of determination of the 
shower direction spherical angles . The difficulty results from the connection (|32p . since the 
"horizontal" components are estimated quantities with random errors. 

Let us employ again the method of error propagation [13] and calculate the matrix of first 
derivatives of the angles ()28p . This time only two "horizontal" components are independent 
variables. 

It turns out that the error of the zenith angle is strongly rising (as a"^ = Dh/cos^{6)] 
see Figl3|) as the shower direction tends to the horizon. This singularity is the strict consequence 
of relation ([32]): if the ort projection on the horizontal plane approaches unity, the reconstructed 
"vertical" component estimation becomes worse evaluated. 

Furthermore, sometimes the "vertical" ort component, evaluated through the nonlinear re- 
lation (|32p with the "horizontal" components containing some errors, becomes imaginary. Really, 
the estimation of ort projection length onto the horizontal plane can prove to be greater then 



unity due to the errors of those components. It is quite unclear how one should interpret such 
result. 

No such accident takes place using the volumetric EAS goniometer, as all ort components 
are computed on the base of the measured delay periods (all being real numbers) by means of 
the real linear transformation (jllh . It cannot have a complex result. Any error can only affect 
the unit length condition of the spatial ort n. The linear distortion of this type never changes 
its geometrical sense of a direction vector, and both angles can be calculated quite intelligently, 
though possibly with big errors. 

5 Examples 



Figure 2: Three examples of goniometer arranging: 

1) high goniometer - left; 2) low goniometer - center; 3) flat goniometer - right. 

Let us consider three simple examples of EAS goniometers for the purpose of illustration, 
(see Figl2]) The detector number in every example installation is A'^o = 8 with the distances 
between them approximately about 20 m, as in example in reference [9|. If we locate them 
in the apexes of squares with edges of 20 m length (the radii of both tiers are of 14.1m in 
size) and estimate the delay periods' standard error by a rough value of cr = 0.5 m ~ 1.7 ns, the 
aforementioned expression (i25]l gives us Dh = 0.0003125. 

However, let us locate the tiers: 

1. with the distance between them equal to 20 m (the so-called high goniometer). The dis- 
persion of the "vertical" ort component = Dh = 0.0003125; 

2. with the distance between the tiers equal to 4 m (low goniometer ) and the same con- 
figuration of tiers as before. The dispersion of the "vertical" ort component increases: 
D, = 0.0078125; 

3. with all detectors placed in common horizontal plane (flat goniometer ) at the apexes of a 
regular octagon with the same circumcircle radius. As before we obtain Dh = 0.0003125. 
(The "vertical" component is indefinite.) 

The dependencies of angles error estimations on zenith angle value are shown on Figl3]for 
all three cases. As one can see, the standard deviation of zenith angle estimation never exceeds 




Figure 3: The errors dependencies 
of spherical angles estimations on the 
true zenith angle value for all consid- 
ered cases of the EAS goniometer ex- 
amples. 
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the value of 1° for high goniometer; while for low goniometer it somewhat exceeds 10° limit only 
for almost horizontal showers. 

However, the LA AS group [9] (together with the most part of investigators all over the 
world) uses the flat goniometer installation. All detectors are placed on the same horizontal 
plane preventing the possibility of linear estimation of vertical component of EAS direction ort. 

On Figl3] the angle dependencies of error estimations for last flat goniometer example are 
shown, too. The azimuth angle is measured with the previous accuracy as the detectors' number 
and their positions radius have not changed. 

The zenith angle estimation error has grown badly. Practically every close-to-horizon angles 
cannot be measured as there estimations coincide with the horizon within the value of standard 



Figure 4: The probability of 
complex zenith angle estima- 
tion. 
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deviation. The angles in the shaded area on FiglS] correspond to this condition. The flat EAS 
goniometer of the last example is not sensitive for zenith angles larger then 68° . 
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Figure 5: Zenith angle esti- 
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The anticipated probability of dummy complex estimation of zenith angle is shown on 
FiglDfor the flat EAS goniometer considered here, as a function of the true zenith angle value. 
Actually, the estimation can become complex in (nearly) the same shaded area on FigJSj 

The complexity of the zenith angle estimation indicated, caused by the fluctuations of the 
horizontal components' estimators, results in the complexity of the expectation value of the ver- 
tical component estimator for any value of true zenith angle. The events resulting in complex 
estimation of the vertical component of the EAS directional ort are plainly rejected in practice. 
This means the use of the real part of complete estimator ([32]) for the vertical component esti- 
mation. This (real) estimator approaches stochastically the real part of the expectation value of 
complete estimator 5ReM(n2), which does not coincide with the true value of the vertical compo- 
nent. Hence, this estimator on the trimmed sample is biased and inconsistent one. Nevertheless, 
the use of this value within the interval of angles with negligible probability of complex estima- 
tions (FigS]) is defensible as the angle's bias value does not exceed the bounds of one standard 
deviation of the resulting zenith angle (Figl5]). 



6 Conclusion. 

The flat kind of goniometer installation has a number of unpleasant properties. When the 
EAS arrival direction lies far from the zenith, the possibility arises of a complex estimation of 
the zenith angle with no clear interpretation. The standard error of the last angle grows rapidly 
to infinity for EAS arrival directions near the horizon. This behavior results in the assertion 
of an insistent desirability to only use volumetric EAS goniometers, especially for EAS network 



stations with big angular distances between them. Even a smah vertical displacement of part of 
the detectors in a flat goniometer array (i.e. conversion to low EAS goniometer) fundamentally 
changes the angles computation conditions: in no case does any complex result arise and the 
error in zenith angle proves to have a superior limit even for horizontal EAS directions. Finally, 
the proper detectors arrangement provides the isotropy of errors. 
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